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The first-principle approach to the dense state of QCD matter, i.e. the lattice-QCD simulation 
at finite baryon density, is not under theoretical control for the moment. The effective model study 
based on QCD symmetries is a practical alternative. However the model parameters that are fixed by 
hadronic properties in the vacuum may have unknown dependence on the baryon chemical potential. 
We propose a new prescription to constrain the effective model parameters by the matching condition 
with the thermal Statistical Model. In the transitional region where thermal quantities blow up in 
the Statistical Model, deconfined quarks and gluons should smoothly take over the relevant degrees 
of freedom from hadrons and resonances. We use the Polyakov-loop coupled Nambu-Jona-Lasinio 
(PNJL) model as an effective description in the quark side and show how the matching condition is 
satisfied by a simple ansatz on the Polyakov loop potential. Our results favor a phase diagram with 
the chiral phase transition located at slightly higher temperature than deconfinement which stays 
close to the chemical freeze-out points. 
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INTRODUCTION 



XT) 

Exploration of the QCD (Quantum Chromo dynamics) phase diagram, particularly toward a higher baryon-density 
regime, is of increasing importance in both theoretical and experimental sides [IJ. From the theoretical point of 
view, so far, only the lattice-QCD simulation [Tj-Q is the first-principle calculation of QCD at work to explore the 
phase transitions associated with chiral restoration and quark deconfinement. The Polyakov loop $ and the chiral 
condensate (^ip) are the (approximate) order parameters for quark deconfinement and chiral restoration, respectively, 
Qh| which are gauge invariant and measurable on the lattice (though both require renormalization corrections). The 
lattice-QCD simulation is, however, of limited practical use and it works only when the baryon chemical potential /ie 
is sufficiently smaller than the temperature T. For \i^JT > 1 the notorious sign problem prevents us from extracting 
any reliable information from the lattice-QCD data p], [5j. 

The effective model study is an alternative and pragmatic approach toward the phase diagram of dense QCD. Some 
may complain that the model study relies on not QCD directly but on just a model. Results from the model analysis 
are, nevertheless, what we can get at best for the moment. Even within the framework of the model study there are 
several different attitudes. One way for theorists to go is simplify QCD so that it can be solvable without introducing 
t^-j- ■ further approximations. QCD-like models in lower dimensions (such as the 't Hooft model) Q, the strong-coupling 
expansion in the lattice formulation [7j, and the large- N c limit of QCD [8] are typical examples in this direction. 
Here, we shall take another way to proceed into the phase structure. The idea is the following (as schematically 
■ illustrated in Fig. []]); 

1. Construct a model that works for infinitely heavy quarks (m q — >• oo) in such a way that the model respects the 
global symmetry (center symmetry) in the finite-T pure gluonic sector. 

2. Choose a chiral model based on the global symmetry (chiral symmetry) for massless quarks (m q — »• 0) in such 
^ ■ a way that the spontaneous breakdown of chiral symmetry is correctly described. 

3. Interpolate a finite-m g model between above-mentioned two. It is minimally required that the infinite m q limit 
and the vanishing m q limit should recover the above models respectively. 

4. Check if the interpolation is properly chosen or not by comparing the model outputs to available lattice-QCD 
and/or phenomenological data. 

Along this line the model is not necessarily solvable and usually needs some additional approximations. Nevertheless, 
if the item 4 above is taken into account very carefully, one may claim that one is dealing with a phase diagram of 
QCD, not of QCD-like models, in a sense that the situation one is handling is not (1 + 1) dimensions, not g 2 — >• oo, 
and not 7V C — >> oo. Sometimes, to this aim toward the QCD phase diagram, one has to face "dirty" businesses; it is 
often the case that the phase structure might be significantly changed by uncontrollable model parameters, which one 
can take in twofold ways — pessimists would be disappointed and say that the model cannot predict anything, and 
on the other hand, optimists would be delighted and say that the model has clarified a non-trivial role played by the 
model parameter in understanding the phase diagram. 

Let us briefly explain how to implement the items 1 and 2. In the absence of particles transforming in the color 
fundamental representation, the genuine gauge symmetry possessed by the pure gluonic theory is S\J(N C )/Z(N C ). If 



2 



consistency? 
Data 





Extrapolation 




Model Pure-gluonic models simulated by 

Monte-Carlo lattice calculations 



Interpolation 

This has been done 

in PNJL and PQM Models 



Chiral models satisfying 
low-energy theorems 



n 



00 



m 



phys 



QCD 



FIG. 1: Schematic picture to show how the effective model is constructed as an interpolation between the pure-gluonic theory 



(at m q <C Aqcd and T ~ Aqcd) and the chiral models (at 



w^phys and T < f^). The prediction from the interpolated 



effective model is done as an extrapolation toward some new axis such as the baryon chemical potential, whose consistency 
with available data must be eventually checked. 



one performs the Z(7V C ) transformation on the gauge links, the fields are shifted typically by 27r/(N c a) where a is 
a characteristic scale (lattice spacing). The perturbation theory breaks this Z(iV c ) symmetry but this is practically 
no problem because the shift goes infinity as a — » 0. Furthermore one can generalize the similar procedure onto not 
the individual gauge link but a product of the N T gauge links along a finite extent. Then, the fields are shifted by 
2tt I (N c N T a) = 2irT/N c , which remains sensible in the a — >• limit. In this way the Polyakov loop matrix L is defined, 
that is, L = Yl T U T and the Z(7V C ) symmetry with respect to L is called "center symmetry" which breaks in the 
perturbation theory [9] . 

The expectation value of the traced Polyakov loop, = (trL), is the order parameter for the quark deconfinement 
phase transition in the pure gluonic system. The most intuitive way to understand this comes from the property 
that is related to the free energy gain of a static single quark placed in a hot gluonic medium as $ = exp[—f q /T]. 
Therefore $ = implies f q = oo, meaning that quarks never show up (confinement). Once becomes non- vanishing, 
f q should take a finite value and thermal quark excitations are permitted. The effective action for L or $ has been 
computed perturbatively [lOj but in order to discuss the phase transition from center symmetric to center broken 
phases, one needs a non-perturbative evaluation of the effective action. 

Concerning the chiral dynamics, the model choice could be anything as long as it can correctly describe the 
dynamical chiral symmetry breaking pattern. Then, the chiral properties are almost automatically derived from 
the so-called low-energy (soft-pion) theorems. Of course some details of the phase transition such as the critical 
temperature and the thermodynamic quantities depend on the choice of the chiral model. Because we are interested 
in the phase transition associated with restoration of chiral symmetry, the non-linear representation is inappropriate 
which is based on the symmetry breaking. 

The order parameter for the chiral phase transition is given by the chiral condensate (ipip}. This is simple to 
understand — the quark mass m q and an operator ^ are conjugate to each other, so m q is a source to generate (ipip) 
and in turn is a source to generate the dynamical mass that breaks chiral symmetry. There are well-established 
chiral models such as the Nambu-Jona-Lasinio (NJL) model and the quark-meson (QM) model. 

The interpolation at the item 3 is the main problem. There is no theoretical justification at all for the existence 
of reasonable interpolation. We can judge how good or how bad it is only through the comparison at the item 4. 
At this point it is already obvious that the Polyakov loop model is not sufficient to access the realistic QCD phase 
transition, though it may capture interesting phenomenological consequences In a similar sense conventional 

chiral models are not good enough to draw the QCD phase diagram even though they are usually designed to be a 
good description of hadronic properties in the vacuum [l2|. To address the QCD phase transitions the first test for 
the validity of the model description should be the consistency check with the known properties available from the 
lattice-QCD simulation at T ^ and /ie = 0; the finite-T behavior of two order parameters, and (V^)? an d the 
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FIG. 2: Chemical freeze-out points taken from Refs. [T3, [Tj|. The red and blue (upper and lower) bands represent the regions 
where the entropy density s and the baryon number density n, respectively, increase quickly from 0.3 to 0.8 in the unit of free 
quark-gluon values, Sf ree and nf ree (see Eq. QJ). 



thermodynamic quantities such as the pressure, the internal energy density, the entropy densit y, e tc. 

Along this line the Polyakov-loop coupled chiral models such as the PNJL (Polyakov-NJL) jl3l . [lij and the PQM 
(Polyakov-QM) [TBI, [l6| models are quite successful to treat both order parameters on the equal footing. Besides, the 
Polyakov loop potential U[<f>] is determined by the lattice data in the pure gluonic theory, namely, by the Polyakov 
loop 3>(T) and the pressure p(T) as functions of T. This means that the PNJL and PQM models include the pressure 
contribution from gluons as well as quarks, so that the models are able to deal with the full thermodynamics which 
are to be compared with the full lattice-QCD simulation. The important point is that the dynamics of transverse 
gluons Af is also under the control of the deconfinement order parameter <E> and thus is to be encompassed in the 
parametrization of the Polyakov loop potential while the Polyakov loop itself is expressed in terms of the 

longitudinal gluon A4. 

Here we would like to emphasize that the success of the PNJL and PQM models is far beyond the fitting physics. 
Model parameters are fixed separately in two regions, i.e. in the pure-gluonic theory (at m q <C Aqcd and T ~ Aqcd) 
and in the chiral models (at m q < m p h ys and T < f n ). The interpolation procedure does not involve any further 
fitting. It is a highly non-trivial discovery that there exists a reasonable way to make an interpolation fairly consistent 
with the full lattice-QCD data. 

The next step one should think of is the prediction from the model. This is done by an extrapolation of the model 
toward some new axis such as the baryon chemical potential. By now there are different versions of the "QCD phase 
diagram" drawn in this way by means of the PNJL and PQM models with different parameter tunings [l3l4l7j. If we 
go into small details, there are many places where we can talk about the difference. Here we shall limit ourselves to 
look at the difference mainly in the behavior of the Polyakov loop or the deconfinement (crossover) transition line. 
Some of the model results may be close to the true answer, and some may not. We must have a guiding principle to 
select out which is preferred and which is not. The available and reliable data at finite baryon density is, however, 
extremely limited. In what follows we shall elucidate the idea and find that the naive extrapolation from the PNJL 
and PQM models is not acceptable. To see this, we will explain the results from the thermal Statistical Model in the 
next section. 
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II. THERMODYNAMICS FROM THE STATISTICAL MODEL 



Regarding the QCD phase diagram at finite T and /ib useful information is quite limited. The lattice-QCD at finite 
density is being improved, but still different techniques to circumvent the sign problem lead to different results. Only 
the chemical freeze-out points in the heavy-ion collisions are experimental hints about the phase diagram. Although 
the freeze-out points shape an intriguing curve on the fiB~T plane, as plotted by error-bar dots in Fig. [2j one should 
carefully interpret it. 

The freeze-out points are not the raw experimental data but an interpretation through the thermal Statistical Model 
[HI, [HI . In this model the grand canonical partition function is given by contributions from the non-interacting gas of 
hadrons and resonances. In view of the fact that the Statistical Model is such successful to fit various particle ratios 
with fiB and T only (/ig, /i s , and fi c are determined by the collision condition), it should be legitimate to take the 
freeze-out points for experimental data, which in turn validates the Statistical Model (though why it works lacks for 
an explanation from QCD). 

Let us proceed by further accepting that the Statistical Model is a valid description of the state of matter until the 
freeze-out curve or even slightly above. It is then a straightforward application of the Statistical Model to estimate 
thermodynamic quantities such as the pressure p, the entropy density s, the baryon number density n, etc. We here 
utilize the open code THERMUS ver.2.1 to calculate s and n at various T and [1b Hq[. From now on the Statistical 
Model analysis specifically means the use of THERMUS. 

Figure [2] shows the chemical freeze-out points taken from Refs. pH Qil, on which s and n are overlaid. For 
convenience we normalized these quantities by 



These are the entropy density and the baryon number density of free massless N% — 1 gluons and N C N? quarks. 

Here we note that, in drawing Fig. [2j we have intentionally relaxed the neutrality conditions for electric charge 
and heavy flavors and simply set \iq = /i s = /i c = 0. We have done so in order to make it possible to compare the 
results from the Statistical Model to the chiral effective model in later discussions. [We note that one can force the 
chiral model to satisfy neutrality but it would be technically involved [2l|.] Nevertheless, we would emphasize that 
the neutrality conditions have only minor effects on the bulk thermodynamics and make only small differences in any 
case. We should also mention that we used Eq. (pQ) with N c = iVf = 3. The choice of Sf ree and nf ree (and relevant TVf) 
is arbitrary and the following discussions do not rely on this particular choice, for we will use Sf ree and nf ree just as 
common denominators to display the Statistical Model results and the PNJL model results. 

The Statistical Model cannot tell us about the QCD phase transitions. Still, Fig. [2] is already suggestive enough. 
We can clearly see the thermodynamic quantities from the Statistical Model blowing up in a relatively narrow region. 
The red and blue (upper and lower) bands indicate the regions where s/sf ree and n/nf ree , respectively, grow quickly 
from 0.3 to 0.8. In the Hagedorn's picture [22| this rapid and simultaneous rise in s and n has a natural interpretation 
as the Hagedorn limiting temperature above which color degrees of freedom is liberated, i.e. color deconfinement. 

The idea here is to make use of the thermodynamics from the Statistical Model as shown in Fig. [2] to judge if the 
Polyakov-loop coupled chiral models work fine at finite density. We also make an important remark that this idea 
can be effective only up to about [1b < 400 ~ 600 MeV because the chemical freeze-out points start dropping down 
steeply in this density region, which suggests an onset of some new form of matter; an example of such possibilities 
is quarkyonic matter [23|. 



Figure [2] is useful to make a guesstimate about the deconfinement boundary but we can deduce no information 
about the chiral property. This is because the thermal Statistical Model needs no medium modification driven by 
chiral restoration. So, to address the QCD phase transitions and associated boundaries, we must find a way to connect 
the thermodynamics in Fig. [2] to the order parameters $ and (V^)- Here let us go into details of the chiral effective 
model. 

It is crucial to adopt the Polyakov-loop coupled model here because the entropy density should contain contributions 
from gluons which are taken care of by the Polyakov loop potential U[<&\. The PNJL model that we use below is 
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FIG. 3: Entropy density normalized by Sf ree (from 0.3 to 0.8) in the Statistical Model (bottom band with red color; same as 
shown in Fig. [2} and that in the PNJL model with a choice To = 200 MeV (top band with green color). The blue band between 
two represents the results with the ansatz (0. 



defined with the following potential; 

w[M] = t 4J a(T) < 



+ b(T) In 1 - 6$$ + 4($ 3 + $ 3 ) - 3($$) 5 



(2) 



with a(T) = a + ai(T /T) + a 2 (T /T) 2 and b(T) = b 3 (T /T) 3 . There are thus five parameters one out of which 
is constrained by the Stefan-Boltzmann limit, i.e. U — » — (87r 2 /45)T 4 at $ = $ = 1 in the T — >• oo limit. These 
parameters are fixed by the pure-gluonic lattice data as ao = 3.51, a\ = —2.47, a 2 = 15.2 , fog = —1.75 [Hj], and 
To = 270 MeV from the deconfinement temperature of first order in the pure-gluonic theory [24]. It is important to 
note that only To is a dimensional parameter, so that the energy scale is set by this To alone. 

In addition the NJL sector of the PNJL model has five more parameters in the three-flavor case [13| appearing in 
the mean-field thermodynamic potential; 



r A d 3 p 



^njl = 9s{{uu) 2 + (dd) 2 + (ss) 2 ) + 4g D (uu)(dd)(ss) - 27V C ^ J e t (p) 



2T 



v f d3p f 

Z?J (27T)3\ 



In det 



I _j_ £ e -(ei(p)-v><d/T 



In det 



1 + T t e" (£l(p)+/i ' 



(3) 



where the energy dispersion relation £i(p) depends on the flavor index i as Si(p) 2 = p 2 + M 2 and the constituent quark 
masses are M u = m u — 2gs(uu) — 2g]j(dd}(ss} and so on. The model parameters are then the light and heavy quark 
masses m u d and ra s , the momentum cutoff A, the four-fermionic interaction strength g$, and the U(l)A-breaking 
six-fermionic interaction strength g^, which are all fixed by the pion mass m^, the kaon mass m^, the eta-prime 
mass rrirjf, the pion decay constant f^, and the chiral condensate (ipip} [12| . 

In the presence of dynamical quarks, if we keep using the pure-gluonic critical temperature To = 270 MeV, the 
simultaneous crossover temperature of deconfinement and chiral restoration becomes over 200 MeV, which is too high 
as compared to the lattice-QCD value. It is nicely argued in Ref. [15] that the back-reaction from quark loops affects 
the mass scale % which changes from T = 270 MeV for N { = down to T = 208 MeV for N { = 2 and T = 187 MeV 
for Nf = 2 + 1 1151. Here we choose to use To = 200 MeV for calculations at /is = throughout. 

In Fig. [3] we show the entropy density calculated in the mean- field PNJL model with To = 200 MeV fixed, in the 
same way as in the Statistical Model drawn in Fig. [2j The bottom (top) band in red (green) color is the result from 
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the Statistical Model (PNJL model). From the figure it is obvious that the PNJL model is not consistent with the 
Statistical Model even at the qualitative level. With the properly scaled To from 270 MeV down to 200 MeV, the 
blow-up behavior in s from the Statistical Model can be smoothly connected to the PNJL model description only 
in the region up to /iB < 400 MeV. The curvature of the band as a function of /is is significantly different; the 
PNJL model result is too flat horizontally and the green band eventually takes apart from the red region where the 
Statistical Model breaks down. 



Such a manifest discrepancy from the Statistical Model is a critical drawback of the PNJL model. The situation 
is not changed even in the PQM model as long as To is a constant. To make the entropy density at /iB ^ 400 MeV 
get saturated earlier as is the case in the Statistical Model, quark degrees of freedom must be released at smaller 
temperature than predicted by the PNJL model. 

One can imagine how this drawback occurs in the PNJL model study; the energy scale in the Polyakov loop potential 
is specified by the dimensional parameter To that may differ depending on the effects of T and /iB in the quark sector. 
We have shifted To from 270 MeV down to 200 MeV, through which we have incorporated the scale change induced 
by TVf quarks at finite T. In this way we may well consider that Tq should decrease with increasing /iB as pointed out 



Our idea proposed here is to make use of Fig. [3] to fix To(/iB) for consistency with phenomenology. One can pick 
up other thermodynamic quantities than the entropy density like the internal energy density, which would anyway 
make little change in the final result. Besides, the choice of the entropy density is most natural because it counts the 
effective degrees of freedom and thus is a sensitive quantity to probe deconfinement. In Ref. [18[ the freeze-out curve 
is parametrized as 



with the fitting result a = 166(2) MeV, b = 1.39(16) x 10~ 4 MeV -1 , and c = 5.3(21) x 10" 11 MeV -3 . Because the 
behavior of the entropy density must be dominantly controlled by deconfinement, we postulate that the change in T 
is to be correlated with Tf(/iB). We see that the freeze-out points and the entropy band in Fig. [2] have roughly same 
curvature indeed. Let us simply use same b and make an ansatz as 



which yields the blue band in the middle of Fig. [3l [To prevent unphysical negative To for large /iB we set a threshold 
at 10 MeV so that Tq > 10 MeV. Hence, the results at T < 10 MeV are not meaningful.] We see at a glance that 
the results from this modified PNJL model have a reasonable overlap with the Statistical Model results in the whole 
density region as plotted. 

At this point one might have thought of several questions. First, the ansatz ([5]) might look ad hoc, but we 
point out that our choice happens to be very close to the independent argument in Ref. [15], in which the /iB- 
dependence has been estimated from the running coupling constant as To(/iB) = T r e — 1// ^ Q:o ^^ B ^ which is expanded 
to be To(/iB)/To c± 1 — 2.1 x 10 _2 (/iB/To) 2 + In the perturbative manner one can also understand how the 
/iB-dependence enters the Polyakov loop potential which consists of the closed loop of dressed gluon propagator. The 
quark-anti-quark polarization diagrams inserted in the gluon propagator generate the back-reaction dependent on 
/iB- There is another phenomenological ansatz for the /iB-dependent U[<f>] [26|. Second, one might wonder if the 
energy scale in the quark (NJL) sector should be modified as well. Such modification is not necessary, however. This 
is because, as we have mentioned, the Statistical Model requires no modification associated with the chiral dynamics, 
which strongly implies that we do not have to introduce /is dependent changes in the NJL parameters. Third, the 
ansatz (|5j) has terms only up to the quadratic order. This means that we cannot go to high-density regions with 
/iB C To. This is indeed so and we have actually truncated higher-order terms in Eq. ((5]). In any case, as we have 
noted before, the idea of the entropy matching holds only up to /iB < 400 ~ 600 MeV, and we should not take the 
results in the high-density region seriously. 



IV. MATCHING PRESCRIPTION 



in Ref. [15|. 



T f (/i B ) = a - bfjb% - c/x| 



(4) 




(5) 



V. 



PHASE DIAGRAM 



Now we get ready to proceed to the possible QCD phase diagram that is fully consistent with the Statistical Model 
thermodynamics in Fig. [2j Using the standard computational procedure of the mean-field PNJL model we can solve 
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FIG. 4: A figure taken from Ref. [2?]]. Phase boundaries associated with deconfinement (blue band) and chiral restoration (red 
band). Each band represents a region where the (normalized) order parameter develops from 0.4 to 0.6. 

<!> and (ipip) as functions of T and /iB, from which the phase boundaries of deconfinement and chiral restoration are 
located. 

Figured! (the central result of Ref. [27j) shows the phase diagram from the modified PNJL model. The blue (red) 
band is a region where the Polyakov loop (normalized light-quark chiral condensate) increases from 0.4 to 0.6. In 
contrast to the old PNJL model, the new results indicate that the chiral phase transition is almost parallel to and 
entirely above the deconfinement, which agrees with the situation considered phenomenologically in Ref. [25j. We 
have found the critical point [H, [29| at (/ie,^) — (45 MeV, 330 MeV), but should not take it seriously since its 
location is beyond the validity region of the current prescription. 

VI. CONCLUSIONS 

It is an intriguing observation that the chiral phase transition occurs later than deconfinement. This is quite 
consistent with the Statistical Model assumption. In the Statistical Model the hadron masses are just the vacuum 
values and any hadron mass/width modifications are not considered, which would be a reasonable treatment only if 
the chiral phase transition is separated above the Hagedorn temperature. Under such a phase structure, besides, our 
assumption of neglecting /ie-dependence in the NJL-model parameters turns out to be acceptable in a similar sense 
as the Statistical Model treatment. This can be understood from the fact that the NJL part yields the hadron masses 
in the vacuum which are intact in the Statistical Model. 

The failure of the standard PNJL model is attributed to baryonic degrees of freedom which is missing; the singlet- 
part of the thermal excitation in the PNJL model can be translated into an expression in terms of baryons as 

f d 3 k Nc{ ^ k 2 +M 2_^ )/T 1 f d 3 k f ( yV 2 +M&- MB )/T / fi x 
J (2^) 3 J (2tt)3 ' W 

where Mn = N c M q: {1b = N c {i q , and k' = N c k with which a factor l/N 3 appears from the integration measure. 
Therefore, the PNJL model significantly underestimates the baryonic excitations by 1/N%. Hence, one may say 
that a modification made in U[$>] by hand stems, in principle, from confinement effects, which can be presumably 
parametrized by the Polyakov loop. This idea is reminiscent of the treatment of transverse gluons. It is an important 
question how our phenomenological input (j5j) is validated/invalidated from the first-principle QCD calculation. This 
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may be answered by future developments in the functional renormalization group method [3C 
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